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Abstract 

We consider the quantum-mechanical problem of a relativistic Dirac particle mov- 
ing in the Coulomb field of a point charge Ze. In the literature, it is often declared 
that a quantum-mechanical description of such a system does not exist for charge val- 
ues exceeding the so-called critical charge with Z = a^ 1 = 137 based on the fact 
that the standard expression for the lower bound state energy yields complex values 
at overcritical charges. We show that from the mathematical standpoint, there is no 
problem in defining a self-adjoint Hamiltonian for any value of charge. What is more, 
the transition through the critical charge does not lead to any qualitative changes in 
the mathematical description of the system. A specific feature of overcritical charges 
is a non uniqueness of the self-adjoint Hamiltonian, but this non uniqueness is also 
characteristic for charge values less than the critical one (and larger than the subcrit- 
ical charge with Z = (y / 3/2)a _1 = 118). We present the spectra and (generalized) 
eigenfunctions for all self-adjoint Hamiltonians. The methods used are the methods 
of the theory of self-adjoint extensions of symmetric operators and the Krein method 
of guiding functionals. The relation of the constructed one-particle quantum mechan- 
ics to the real physics of electrons in superstrong Coulomb fields where multiparticle 
effects may be of crucial importance is an open question. 

Keywords: Dirac Hamiltonian, Coulomb field, self-adjoint extensions, spectral analysis. 



1 Introduction 

It is common knowledge that the complete sets of solutions of relativistic wave equations 
(like Klein-Gordon equation, Dirac equation, and so on) when used in quantizing the corre- 
sponding (scalar, spinor, and so on) free fields allow interpreting the corresponding quantum 
theories in terms of particles and antiparticles pQ. The space of quantum states of each 
such a field is decomposed into sectors with a definite number of particles (the vacuum, 
one-particle sector, and so on), each sector is stable under time evolution. A description of 
the one-particle sector of a free QFT can be formulated as a relativistic quantum mechanics 
where the corresponding relativistic wave equations play the role of the Schrodinger equation 
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and their solutions are interpreted as wave functions of particles and antiparticles. In QED 
(and some other models), the concept of the external electromagnetic field is widely and 
fruitfully used. It can be considered an approximation in which a "very intensive" part of 
the electromagnetic field is treated classically and is not subjected to any back reaction of 
the rest of the system. The Dirac equation with such a field plays an important role in QED 
with an external field. Of special interest are the cases where an external field allows exactly 
solving the Dirac equation. There are a few of such exactly solvable cases of physically 
interesting external electromagnetic fields, see, e.g., p]. They can be classified into groups 
such that the Dirac equations with fields of each group have a similar interpretation. 

The constant uniform magnetic field, the plane-wave field, and their parallel combination 
form a first group, the fields of this group do not violate the vacuum stability (do not create 
particles from the vacuum). The exact solutions of the Dirac equation with such fields 
form complete systems and can be used in the quantization procedure providing a particle 
interpretation for a quantum spinor field in the corresponding external background. This 
allows constructing an approximation where the interaction with the external field is taken 
into account exactly, while the interaction with the quantized electromagnetic field is treated 
perturbatively. Such an approach to QED with external fields of the first group is known 
as the Furry picture, see, e.g., 0, OQ. In the Furry picture, the state space of the quantum 
theory of the spinor field with the external fields is decomposed into sectors with a definite 
number of particles, each sector is stable under time evolution similarly to the zero external- 
field case. The description of the one-particle sector also can be formulated as a relativistic 
quantum mechanics [4j. We note that the solutions of the Dirac equation with an uniform 
magnetic field provide a basis for the quantum synchrotron radiation theory [5], and the 
solutions of the Dirac equation with the plane-wave field are widely used for calculating the 
quantum effects when electrons and other particles of spin one-half move in laser beams [6]. 

A uniform electric field and some other electromagnetic fields violate the vacuum stability. 
A literal application of the above approach to constructing the Furry picture in QED with 
such fields is fails. However, it was demonstrated that existing exact solutions of the Dirac 
equation with electric-type fields form complete sets and can be used for describing a vsriety 
of quantum effects in such fields, in particular, the electron-positron pair production from 
the vacuum [TJ. Moreover, these sets of solutions form a basis for constructing a generalized 
Furry picture in QED with external fields violating the vacuum stability, see [8]. It should be 
noted that the one-particle sector in such external fields is unstable under time evolution, and 
therefore, the corresponding quantum mechanics of a spinning particle cannot be constructed 
in principle. 

A study of the Dirac equation with a singular external Aharonov-Bohm field, and with 
some additional fields revealed problems of a new type. Although some sets of exact solutions 
of the Dirac equation with such fields can be found, the problem of the completeness of these 
sets arises. This problem is related to the fact that the Dirac Hamiltonian with the singular 
external Aharonov-Bohm field should be additionally specified for it can be treated as a self- 
adjoint (s.a. in what follows) quantum-mechanical operator. It can be shown (for a review, 
see [9]) that in this case, there exists a family of s.a. Hamiltonians which are constructed 
by methods of the theory of s.a. extensions of symmetric operators dating back to von 
Neumann [TO]. Each s.a. Hamiltonian yields a complete set of solutions which can be used 
for constructing the Furry picture in QED with the singular external Aharonov-Bohm field 
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(this field does not violate the vacuum stability). 

The Dirac equation with the Coulomb field, and with some additional fields, has always 
been of particular interest. The Coulomb field is even referred to as a "microscopic external 
field" to underline its qualitative distinction from the above-mentioned external fields which 
are sometimes referred to as "macroscopic" ones. Until recently, the commonly accepted 
view on the situation in the theory was the following. The Dirac equation for an electron 
of charge^ — e in an external Coulomb field created by a positive point-like electric charge 
Ze of a nucleus of atomic numbeJl Z < 1/a = 137 is solved exactly, has a complete 
set of solutions, and allows constructing a relativistic theory of atomic spectra which is in 
agreement with experiment [II] . This field does not violate the vacuum stability, therefore, 
the Furry picture can be constructed, and there exists the relativistic quantum mechanics of 
the spinning particle in such a Coulomb field. As for the Dirac equation with the Coulomb 
field with Z > 137, it was considered inconsistent and physically meaningless [T5] [T3]. One 
of the standard arguments is that the formula for the lower IS1/2 energy level, 



formally gives imaginary eigenvalues for the Dirac Hamiltonian with Z > 137. On the 
one hand, the question of the consistency of the Dirac equation with the Coulomb field 
with Z > 137 has a pure theoretical (mathematical) interest, on the other hand, it is 
concerned with the question of the electron structure of atoms of atomic numbers Z > 1/a 
and especially of atoms with nuclei of supercritical charge Ze > 170e . The latter question is 
of fundamental importance. The formulation of QED cannot be considered really completed 
until an exhaustive answer to this question is given. Although nuclei of so large electric 
charges can hardly be synthesized^], the existing heavy nuclei can imitate the supercritical 
Coulomb fields at collision. Nuclear forces can hold the colliding nuclei together for 10~ 19 s 
or more. This time is enough to reproduce effectively the experimental situation where 
the electron experiences the supercritical Coulomb field [15]. Several groups of researchers 
attacked the problem of the behavior of the electron in the supercritical Coulomb field, 
see [15], [16]. The difficulty of the imaginary spectrum in the case of Z > 137 was attributed 
to an inadmissible singularity of the supercritical Coulomb field for a relativistic electron 
@. It was believed that this difficulty can be eliminated if a nucleus of some finite radius 
R is considered. It was shown that with cutting off the Coulomb potential with Z < 170 
at a radius R ~ 1,2 x lCT 12 cm, the Dirac equation has physically meaningful solutions 
[17] . But even in the presence of the cut off, another difficulty arises at Z ~ 170. Namely, 
the lower bound state energy descends to the upper boundary E = —mc 2 of the lower 
continuum, and it is generally agreed that in such a situation, the problem can no longer 
be considered a one-particle one because of the electron-positron pair production, which 
in particular results in a screening of the Coulomb potential of the nucleus. Probabilities 
of the particle production in the heavy-ion collisions were calculated in the framework of 

1 e = 4, 803 x 1Q~ 10 CGSE is the magnitude of the electron charge. 
2 a = e 2 /he is the fine structure constant. 
3 At present, the maximum Z = 118. 

4 An equation for the radial components of wave functions has the form of the nonrelativistic Schrodinger 
equation with an effective potential with the r~ 2 singularity at the origin, which is associated with "a fall 
to the center" . 
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this conception [15J . Unfortunately, experimental conditions for verifying the corresponding 
predictions are unavailable at present. 

In this paper, we turn back to the problem of the consistency of the Dirac equation 
with the Coulomb field with no cut off and with arbitrary nucleus charge values( with 
arbitrary Z) . Our point of view is that the above-mentioned difficulties with the spectrum for 
Z > 137 do not arise if the the Dirac Hamiltonian is correctly defined operator. We 

present a rigorous treatment of all the aspects of this problem including a complete spectral 
analysis of the model based on the theory of s.a. extensions of symmetric operators and the 
Krein method of guiding functionals. We show that from the mathematical standpoint, the 
definition of the Dirac Hamiltonian cLS cL S.cL. operator for arbitrary Z presents no problem. 
What's more, the transition from the noncritical charge region to the critical one does not 
lead to qualitative changes in the mathematical description of the system. A specific feature 
of the overcritical charges is a non uniqueness of the s.a. Dirac Hamiltonian, but this non 
uniqueness is characteristic even for charge values less than the critical one (and larger than 
the subcritical value with Z = (v / 3/2)a _1 = 118). Presenting a rigorous treatment of 
the problem, we also compare it with the conventional physical approach to constructing a 
quantum-mechanical description of the relativistic Coulomb system. It turns out that many 
of mathematical results can be obtained from physical considerations except the important 
property of completeness for the eigenf unctions of the Hamiltonian. The obtained complete 
sets of solutions can be used to construct the Furry picture in QED. However, it is unclear 
whether the neglect of the radiative interaction in such a Furry picture is a good zero 
approximation to describing quantum effects in QED with the Coulomb field with no cut 
off and with arbitrary nucleus charge values. In other words, a relevance of the constructed 
quantum-mechanics with an energy spectrum unbounded from below to the real physics 
of an electron in supercritical Coulomb fields where multiparticle effects may be of crucial 
importance is an open question. 

The paper is organized as follows. In Sec. [2J we present basic facts and formulas clar- 
ifying the formulation of the problem and reduce the problem of constructing a s.a. Dirac 
Hamiltonian with an external Coulomb field in the whole Hilbert space to the problem of 
constructing s.a. one- dimensional radial Hamiltonians. In Sec. [3j we cite expressions for 
the general solution of the radial equations and some particular solutions of these equations 
used in the following. In Subsec. 14.11 and Appendix 14.5.21 we outline procedures for con- 
structing s.a. extensions of symmetric differential operators and their spectral analysis. In 
Subsecs. I4.2H4.51 we construct s.a. Dirac Hamiltonians in all four possible charge regions 
and find their spectra and the corresponding complete sets of eigenfunctions. 

2 Setting the problem 

We consider the Dirac equation for a particle of charge q\ moving in the external Coulomb 
field of a point-like charge qi\ for an electron in a hydrogen-like atom, we have qi — — e , 
q2 = Ze. We choose the electromagnetic potentials for such a field in the form 



A = ^, A k = 0. 
r 
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The Dirac equation with this field, being written in the form of the Schrodinger equation, 

where \I/ (x) = {ip a (x)} is a bispinor and the Dirac Hamiltonian H is given by 

q ( m-q/r crp 



r 



H = cp + m(3-^=[ y . , (1) 

1 crp — m — g/r ' 



m is the fermion mass, p = (p k = —id/A , and q = —q±q2] for an electron in a hydrogen-like 
atom, we have q = Za. For brevity, we call the coupling constant q the charge. We restrict 
ourselves to the case of q > 0, because the results for the case of q < can be obtained by 
the charge conjugation transformation. 

At this initial stage of setting the problem, the Hamiltonian H and other operators are 
considered as formally s.a. differential operators, or s.a. differential expressions^, as we will 
say [TS] , which is denoted by the turned hat V above the corresponding letter. They become 
quantum-mechanical operators after a specification of their domains in the Hilbert space 7i 
of bispinors ^(r), 

4 

H = 2J T~ia, 7~La = L 2 (R 3 ) 

a=l 

then the symbol V is replaced by the conventional symbol A over the same letter. In what 
follows, we distinguish differential expressions / and operators / and call / the operator 
associated with the differential expression /. 

The purposes of this paper are constructing a s.a. Hamiltonian H associated with H, 
which primarily means indicating a domain of H, and then finding its spectrum and eigen- 
functions. 

We note that in the physical literature, the eigenvalue problem is conventionally con- 
sidered directly in terms of the s.a. differential expression H as the eigenvalue problem for 
the differential equation H^e^) = E^/eM, the stationary Dirac equation, without any 
reference to the domain of the Hamiltonian_|. It is solved by separating variables based on 
the rotation symmetry of the problem. The rotation symmetry is conventionally treated in 
terms of s.a. differential expressions as follows. 

The Dirac Hamiltonian H formally commutes with the total angular momentum J = 
L + |S, where L = [r x p] is the orbital angular momentum operator and |S is the spin 
angular momentum operator, and the so-called spin operator K, 

# = 0[l+(Et)] = (* \ ),*=1+(<tL). 



3 We use the bold letters for there-vectors and the standard representation for 7-matrices where 



CT 

<x 



and er = (er 1 , a 2 , a 3 ) are the Pauli matrices. We use the notation crp =a k p k , err =<j k x k , and so on. We set 
h = c = 1 in what follows. 

6 S.a. by Lagrange in the mathematical terminology, or formally s.a. in the physical terminology. 

7 In fact, some natural domain is implicitly implied. 
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The differential expressions H, J 2 , J3, and K are considered a complete set of commuting 
operators, which allows separating the angular and radial variables and reducing the to- 
tal stationary Dirac equation to the radial stationary Dirac equation with a fixed angular 
momentum, its z-axis projection, and a spin operator eigenvalue. 

We here present a treatment of the problem which is proper from the functional analysis 
standpoint. We construct a s.a. Hamiltonian H based on the theory of s.a. extensions of 
symmetric operators and on the rotation symmetry of the problem. This means that we 
first define a rotationally invariant symmetric operator associated with s.a. differential 
expression H (JT]), which is rather simple, and then find its rotationally invariant s.a. exten- 
sions. Because the coefficient functions of H are smooth out of the origin, we choose the 
space of smooth bispinors with a compact support for the domain D H ( ) of H^°\ To avoid 
troubles with the 1/r singularity of the potential at the origin, we additionally require that 
all bispinors in D H (o) vanish near the origin^. The operator H^°hs thus defined by 

/wo) / D H(0) = {* (r) : a (r) E D(R 3 ); ^ Q (r) = 0, r eU £ }, 
■ \ fl-(o)*( r ) = m(r), 

where D(R 3 ) is the space of smooth functions in i? 3 with a compact support and U e is some 
vicinity of the origin which is generally different for different bispinors. The domain D H (o) 
is dense in H, D H { ) = H, and the symmetricity of is easily verified by integrating by 
parts. 

We now take the rotational invariance into account. The operator evidently com- 
mutes with the s.a. angular momentum operator J = {Jk} and the s.a. operator K associ- 
ated with the respective differential expressions J and K. The operators Jk are defined as 
generators of the unitary representation of the rotation group Spin (3) in the Hilbert space 
Ji. The Hilbert space 7i is represented as a direct orthogonal sum, 

n = Y^ n ^ = 1/2,3/2,..., c = ±i, (2) 

of subspaces T~lj£. The subspaces TCj^ reduce the operators J 2 , Jk, and 

j 2 = ^® J 2 f; j k = ^® J fciC) K = Kj£. 

i.C j,C i>C 

In its turn, the subspaces Ttj^ are finite direct sums of the subspaces 

«i,C= M=-j,-j + l,...,j. (3) 

M 

The subspace 7~tj,M,( is the subspace of bispinors *&j,M,c ( r ) °f the form 



We thus avoid troubles associated with a behavior of wave functions at infinity. 
9 Strictly speaking, we thus leave room for <5-like terms in the potential. 

10 This means that the operators J 2 , J3, and K commute with the projectors to the subspaces see 
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where ^^(fl, f) are spherical spinors and / (r) and g (r) are radial functions (the factors 
1/r and % are introduced for convenience). The subspaces Tij,M,c are the eigenspaces of the 
operators J 2 , J 3 , and 

J 2 ^- M ,c (r) = j(j + l)*i,M,c (r) , i 3 %M,c W = M%, M , 
^• M>c (r) = -C(j + l/2)^. M)C (r), 

and evidently reduce the operators J 2 , J 3 , and In the physical language, decomposition 
© corresponds to the expansion of bispinors ^ (r) G in terms of the eigenfunctions of 
the commuting operators J 2 , J 3 , and K, which allows separating variables in the equations 
for eigenfunctions. 

We note that the reductions J\,2jx or " the operators Ji^ to the subspaces are bounded 
operators. 

The following fact is basic for us. Let £ 2 (0, 00) be the Hilbert space of doublets F(r), 

F(r) = ( / < r » 

with the scalar product 

poo /*oo r 

(F 1; F 2 ) = / drF+ (r) F 2 (r) = / rfr 7TW/ 2 (r) + ^fej (r) 
Jo ■/ L 

such that £ 2 (0, 00) = L 2 (0, oo) e L 2 (0, 00). Then formula flU and the relation 

II%m,cII 2 = /dr*W)*W r ) = / W(0| 2 +l<7(0| 2 ] 
show that the Hilbert space Tlj,M,c * s unitary equivalent to the Hilbert space £ 2 (0, 00): 

the explicit form of the unitary operator U is defined by fl4]). 

The rotational invariance of is equivalent to the following statement. 

1) The subspaces H.j,M,c, reduce this operator, such that is represented as a direct 
orthogonal sum of its parts H<$ and that are the reductions of to the respective 

7ij ^ and 7ij m o 

Each part ijj ]^ is a symmetric operator in the Hilbert space Tlj,M,c- Each symmetric 

operator H^^- in the subspace 7^j,m,c evidently induces a symmetric operator ft,^ in the 
Hilbert space £ 2 (0, 00), 

such that ^ = U jM xHj°M,( U i,M,o and ^i°C is given by 

J ' C \ ft<°>F(r) = fc At F(r), 
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where P(0, oo) = D(0, oo) e /}(0, oo), D(Q, oo) is the standard space of smooth functions on 
(0, oo) with a compact support, 

D (0, oo) = {/ (r) : / (r) G C°°, supp/ C [a, b), < a < b < oo} , 

the segment [a, 6] is generally different for different /, and the s.a. radial differential expres- 
sion hj£ is given by 

L c = -io 2 ^- + — a 1 - - + ma 3 , (7) 
dr r r 

where x= C(j + 1/2)0 

2) The differential expression ft,^ and consequently, with taking (jSj) into account, the 
operator ft,^ with fixed j and £ is independent of M. This fact is equivalent to the com- 
mutativity of the operator with the operators J\ and J2, or more precisely, to the 

commutativity of the operator with the operators Ji,2j,(- m the physical terminology, 
hj^ is called the radial Hamiltonian, but strictly speaking, the radial Hamiltonian is a s.a. 
operator hj£ associated with hj^. 

In what follows, by the rotational invariance of any operator / , we mean the fulfilment 
of the following conditions: 

1) the reducibility of this operator by the subspaces TCj,M,c anc ^ therefore, by TCj£, such 
that formula similar to (jSJ) holds for the operator / ; 

2) the commutativity of its parts fj^ with the bounded operators J\,2j,c for fixed j, £. 
Let hj£ be a s.a. extension of the symmetric operator in £ 2 (0, 00). It evidently 

induces s.a. extensions Hj,M£ of the symmetric operators H^ M ^ in the subspaces TC^mx, 

Hj,M,c = U^M^hj^Uj^, (8) 

and the operator Hj^ = M Hj,M,< commutes with Ji,2j,c- Then the closure of the direct 

orthogonal sum 

11 Y. 11 - Y 11 ^ ( ) 

IS cL S.cL. operator in the whole Hilbert space Ti [19], and H is a rotationally invariant 
extension of the rotationally invariant symmetric operator H^°>. 

Conversely, any rotationally invariant s.a. extension H of the initial operator has 
structure ([!]), and the operator hj^ = U^m^H^mx^Zmx in £ 2 (0, 00) is independent of M 
and is a s.a. extension of the symmetric operator h j°j 1^1 

The problem of constructing a rotationally invariant s.a. Hamiltonian H is thus reduced 
to the problem of constructing s.a. radial Hamiltonians hj^. 

In what follows, we operate with fixed j and ( and therefore omit these indices for 
brevity. In fact, we consider the radial differential expressions hj$ as a two-parameter 



We note that X>(0,oo) is dense in £ 2 (0,oo), because, as is known, D (0, 00) is dense in L 2 (0,oo), and 
symmetricity of hj^* is easily ve 
Roughly speaking, this means tt 
M's must be constructed "uniformly 



the symmetricity of hj 3 ^ is easily verified by integrating by parts, which confirms the above assertion. 
12 Roughly speaking, this means that s.a. extensions of the parts H^ M ^ with fixed jand Q and different 
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differential expression h with the parameters q and x (the parameters j and ( enter through 
one parameter x, the parameter m is considered fixed) and similarly treat the associated 
radial operators and h defined in the same Hilbert space £ 2 (0, oo). 



3 General solution of radial equations 

Later on, we need some special solutions of the differential equation 

hF = WF, 

or, which is the same, the system of equations 

f r +?f-(W + m + *)g = 0, 

with an arbitrary complex W; real W are denoted by E and have the conventional sense of 
energy We call Eqs. ffTUj) the radial equations. For completeness, we present the general 
solution of the radial equations following the standard procedure, see, e.g., [13j|14]. We first 
represent f(r) and g(r) as 

f(r) = z T e^ 2 [P (z) + Q(z)} , g(r) = -ihz' 1 % e~^[P{z) - Q{z)} , 

where z = —2iKr , T, A, and K are some complex numbers that are specified below. The 
radial equations then become the equations for the functions P and Q. Setting 

T 2 = x 2 -g 2 , a = r-4^, /3 = 1 + 2T, 

K 

W±m = p±e iLp± , < ip± < 2vr, (11) 



V W + m y p + vt 
we reduce ffTUl) to the system of equations 

F= - x-l mj K) (4 + °)» < 12) 

The first equation in ( Tl2l) is the confluent hypergeometric equation [2D1 El] for Q. 
Let@ T 7^ —n/2, n = 1, 2, then its general solution can be represented as 

Q = A$(a, p; z) + BV(a, /?; z) , (13) 



13 The parameter T is defined by (JTTJ) up to a sign. A specification of T is a matter of convenience. 
In particular, for specific values of charge, we also use a specification of T where T = —n/2, this case is 
considered separately below. 
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where A and B are arbitrary constants, $>(a,/3;z) and ^(a, /3; z) are the known confluent 
hypergeometric functions, for their definition, see [201 I2T] (the function $>(a,/3;z) is not 
defined for (3 = 0, -1, -2, ...)• It follows from eqs. jT2) and (USD that 



P 



a 



x — i(qm/K) 



aW 

A$(a + 1, # z) - B{T + (« + 1; P\ z ) 

K 



The general solution of radial equations (flOj) for any complex W and realm, x, and q is 
finally given by 

/( r ) = 2 T e- z/2 {A[$(a, 0; 2) - a+$(a + 1, # z)} + (a, # 2) + b^{a + 1; # z)}}, 
g(r) = iAz T e' z/2 {A[^(a, (3; z) + a+$(a + 1, /?; 2)] + (a, /?; z) - + 1; /?; 2)]}, 

a± = F : ' h = T? • 

kK — iqm K 

Taking the relation 

$(a + 1, p- -2iKr) = e~ 2iXr $(/3 -a -1,0; 2iKr) 



into account, see [201 El], it is convenient to represent the general solution of radial equations 
( flOl) in the form 



AX(r, T, W) + Bz r e- z/2 [tf (a, /3; 2)$+ - + 1, # ^)^_ 



(14) 



where the doublets &± are 
and the doublet X is 



0- 



X 



[mr] 



$^ = e ^ r $ I T + ±— , 1 + 2T; -2iKr 
iK 



$+(r,T, W) + $-(r,T,W) 



±1 

zA 



m + W 
m-W 



iKr > 



$ T - 



tk 



, 1 + 2T; 2ifsTr 



(15) 



1 



e* Kr $ I T + 



~iK' 



1 + 2T: -2iiTr I - e~ iKr ^ I T 



qW 

Ik' 



1 + 2T; 2iKr 



the doublet w + is one of doublets u± which are used below, 



u± 



1 



We now present some particular solutions of radial equations f llOp which are used in the 
following. 

One of the solutions given by (fT4"j) with A = 1, B = 0, and a specific choice for T is 



U {1) (r;W) = X(r,T,W)\ 



T=T+ ' 



(16) 
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where 

v/x 2 — q 2 , q < |x| 



q 2 — x 2 , q > |x| 



2 



in what follows, we set 7 = yx 2 — g 2 anckr = y q 2 — x 2 . The asymptotic behavior of the 
doublet W) at the origin is given by 

U {1) (r; W) = (mr) T +u + + 0(r T++1 ), r -> 0. (17) 

In the case where T + ^ n/2, n = 1, 2, we also use another solution 

C/ (2) (r;W) = X(r,T,W)| T= _ x+ (18) 

with the asymptotic behavior 

U (2) (r; W) = (mr)~ T +w_ + 0(r~ T++1 ), r -> 0. (19) 

For q 7^ g CJ - = |x| = j + |, i.e. for T + 7^ 0, the solutions U(x)(r;W) and f/(2)( r ;W / ) are 
linearly independent, 

2T 

Wr(f/ (1) ,f/ (2) ) = -±, (20) 

where Wr(F 1; F 2 ) = Fiia 2 F 2 = fig2 — 9if2 is the Wronskian of the doublets F x = ( ^ ] and 

V #1 / 

h 

92 

It follows from the standard representation for the function $ that for real T(T ^ —n/2), 
the functions <3> + and $_ in ( fl5|) are real-entire functions of W, i.e., they are entire in W and 
real for real W = E. It then follows from representations ( fl6l) and ( fl8|) that the respective 
doublets U(i){r; W) and U(2)(r; W) are also real-entire functions of W for real T + = 7. If 
T + is pure imaginary, T + = ia, then £/(i)(r; W 7 ) and {7(2)0"; VT) are entire in W and complex 
conjugate for real W — E, £7(i)(r; £?) = {7( 2 )(r; £7). 

Another useful solution nontrivial for T + ^ n/2, n = 1, 2, is given by ffT4j) with A = 
and a special choice for B: 

V (1) (r; W) = B(W)(mr) r+ e iKr [tf (a, /9; ^)tf+ - (a + 1; /3; z) 0_] , 

_ r(-T + + gwyig) _ 1 

1 ] r(-2T+)(l-a + ) "r(-2T + )B(M0" 1 J 

As any solution, Vm is a special linear combination of {7m and C/( 2 ), 

V(i)(r; = C/ (1) (r; W) + ^o;(W)C/ (2) (r; WQ, (22) 

where 

w(W) = -Wr(t/ (1) ,V( 1) ) = 



2T+r(2T+)r(-T + + qW/iK)(l - a_){2e-™/ 2 K/m)- 2T + _ 



g r(-2T + )r(T + + gW/ilir)(l-a + ) T(-2T. 



(23) 
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We note that if > and r — > oo, the doublet U^(r; W) increases exponentially while 
V(i)(r; W) decreases exponentially (with a polynomial accuracy). 

The doublets Ui 2 ) and Vm are not solutions linearly independent of Un\ at the points 
T + = 7 = n/2 where f/( 2 ) is not defined while Vm vanishes. We need their analogues 
defined at these points and having all the required properties, in particular, real-entirety in 
W. Unfortunately, we can construct such solutions only in some neighborhood of fixed n. 
The corresponding solutions are labelled by the index n. 

According to ( 12T1) . the doublet Vm tends to zero like l/r(— 27) as 7 — > n/2. It then 
follows from (1221) that the doublet U(2) has a singularity of the form T(— 27) at the point 
7 = n/2 and can be represented in a neighborhood of this point as 

C/ (2) (r; W) = r(-2 7 )A n (PV)f/ {1) (r; IV) + C/ n(2) (r; IV), (24) 

where 



g Cb(W) 



7=n/2 



and U n (2){r\ W) has a finite limit as 7 — > n/2 and evidently satisfies radial equations (I10p . 
A direct calculation^ shows that A ra (W) is a polynomial in W with real coefficients, and 
because E/m(r, W) and U^(r, W) are real entire in PV, the doublet £/ n ( 2 )(r; W) is also real 
entire. We thus obtain that the doublet U n ( 2 ) defined by 

U n{2) (r; W) = U m {r- W) - r(-2 7 )A n (VV)C/ (1) (r; W) (26) 

and satisfying the condition 

U n{2 )(r; W) = (mr)~ 7 w_ + 0(r~ 7+1 ), r -> 0, 

is a solution of radial equations which is well-defined in some neighborhood of the point 7 = 
n/2, the point itself included. This solution is linear-independent of £7m, Wr (E/m, E/ n (2)) = 
— 27/g, and is real entire in W. According to relations (I22l - ([2"4"1) . the doublet Vm is repre- 
sented in a neighborhood of the point 7 = n/2 in terms of the finite doublets E/m and U n ^) 
as 

V (1) (r; W) = [1 + ^-u(W)A n (W)]U (1} (r; W) + |- w (W)f/ n(2) (r; W), 

where according to ( 1251) and ( 1231) . the factors 1 + ^-uj(W)A n (W) and u;(W) tend to zero like 
l/r(— 27) as 7 —>■ n/2. This allows introducing the doublet 

K(i)(r; W) = 1 + ± Q{ w )MW) V(i)(r; W) = U (1) (r; W) + ^ n (W)U n(2) (r; W), (27) 

where 

u (W) = U{W) = (2Z) 

n{ ] i + ^mw r(-2 7 )[i + ^(w)A n (iv)]' 1 ] 

which is evidently a solution of the radial equations well-defined in some neighborhood of 
the point 7 = n/2, the point itself included, and exponentially decreasing as r — > 00. The 



14 With the use of the equality T(w + 1) = wr(w) 
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function u n (W) is also well defined in some neighborhood of the point 7 = n/2 and at the 
point itself. We point out that the useful relations 

-K(i)(r; W) = -^—V (1) (r; W), (29) 



U> n (W) " WV ' ' Lo{W) 

' -f r(-2 7 )A n W + tTitT (30) 



hold; strictly speaking, they are meaningful for 7 7^ n/2 and show that the r.h.s.'s are 
continuous in 7 at the points 7 = n/2. 

The doublets U n r 2 )(r;W) and (r; W) are the required analogues of the doublets 
U(2){r; W) and V(i)(r; W) defined in the neighborhood of the point 7 = n/2 and at the point 
itself. 

It remains to consider the special case of q = q C j = j + 1/2, or Y = 0, where the doublets 
£/(i) and Up) coincide while V(i) vanishes. Let C/(i)(r; W 7 ^), [7( 2 )(r; W 7 ^), and V(i)(r;W / | 7 ) 
denote £/( 2 ), and V(i) with 7 7^ 0. Differentiating radial equations ({TO]) for [/(i) with 
respect to 7 at 7 = 0, we can easily verify that the doublet 



dU {l) {r-W) 



<9 7 



7=0 



lim ^(i)(r;^|7)-^(2)(r;^|7) 
7^0 27 



is a solution of these equations with 7 = 0. For two linearly independent solutions of radial 
equations (TT01) with 7 = 0, we choose 



U {1) (r;W) = U {1) (r;W\0), 

U {1) (r;W)=u + + 0{r), r -> 0, (31) 

and 



17$ (r; W0 = ^ ' ' - ^(i)(r; W|0), 

[/^(r; W) = u {0) (r) + O(rlnr), r -> 0, (32) 
where u + and u_ (r) are 

1 "\ ..(0)/_n f ln (™ r ) - 



" + = l c J' "- w n cmW" '■ (33) 



The Wronskian of these solutions is 



Wr([/ (1)) [/$) = -l. (34) 



The both doublets U(i) and uL°) are real entire in W. 

As an analogue of V(\) in the case of 7 = 0, we take the doublet 



V$(r; W) = lim[-r(-2 7 )V( 1) (r; W| 7 )] = (35) 
-^^e iXr [tf (a, 1; -2itfr) + b^(a + 1, 1; -2z^r) a 3 ] ( 
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where 



a 



q CJ W 
%K ' 



Its representation in terms of U(x) and U,J is given by 



W 

m + i(,K' 

(0) 



Qcj- 



(K + im 
~K ' 



r(0) 

(i) 



V$(r-W) 



U$(r; W) + q cj u^\W)U {1) {r- W), 



(2) 



Wr(U$(r-,W),V$(r;W)) 



(0), 



(36) 



1 

9 C j 



ln(2e-^ /2 ir/m) + i/)(-iq cj W/K) + 



2^(1) 



where ip is a symbol of the logarithmic derivative of the T-function. We note that if QW > 

A ' 
a) 



0, the doublet V$(r;W) is square integrable, V,y)(r; W) G £ 2 (0, oo), and exponentially 



decreases as r 



(i) 
oo. 



4 Self-adjoint radial Hamiltonians for different regions 
of charge q 

4.1 Generalities 

In this section, we construct a s.a. radial Hamiltoniarl^l h in the Hilbert space C 2 (0, oo) of 
doublets el. S 9 S.9. extension of symmetric radial operator ([6]) associated with the radial 
differential expression h (JTj) and analyze its spectral properties. 
The extension procedure includes the following steps [22] : 

1) evaluating the adjoint operator h* = (ji^^j and estimating its asymmetricity in terms 

of (asymptotic) boundary values of doublets belonging to the domain of h*. 

2) constructing s.a. extensions h of as s.a. restrictions of the adjoint h* specified by 
some (asymptotic) s.a. boundary conditions at the original. 

It turns out0 that the result crucially depends on the value of the charge q: different 
regions of the charge are assigned different s.a. radial Hamiltonians in the sense that they 
are specified by completely different types of (asymptotic) s.a. boundary conditions. What 
is more, for sufficiently large charg GS, Sb S.9. radial Hamiltonian is defined non-uniquely, such 
that there is a one-parameter family of s.a. Hamiltonians for fixed x and q. Therefore, our 
exposition is naturally divided into subsections related to the corresponding regions of the 
charge; actually, there are four of them. 

For each region, we perform a full spectral analysis of the obtained s.a. Hamiltonians, in 
particular, we find their spectra and (generalized) eigenfunctions. The analysis is based on 
the Krein method of guiding functionals and includes the following steps: 

i) constructing the guiding functional, 

15 We omit indices in the notation of h, h^°\ and h, see the end of Sec. [21 

16 We note that this method of constructing h allows avoiding an evaluation of the deficient subspaces 
and deficiency indices of , the latters are determined by passing. 
17 Actually, this becomes clear at the first step. 
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ii) evaluating the resolvent, 

iii) evaluating the spectral function, 

iv) constructing the so-called inversion formulas that are mathematically rigorous formu- 
las for the Fourier expansion of wave functions with respect to the (generalized) eigenfunc- 
tions of the s.a. Hamiltonian. 

We compare this analysis with heuristic physical considerations based, in particular, on 
the rule of "normalization to 5-function" for eigenf unctions of the continuous spectrum. 

Our first task is constructing a s.a. radial operator h in accordance with the above 
scheme, which mainly reduces to indicating its domain D^, D h @) C Dh Q -D*. It can happen 
that such a domain is non-unique, and it really is for some values of number parameters in 
h. 

The domain for a s.a. differential operator on an interval of the real axis is conventionally 
specified by the so-called s.a. boundary conditions at the ends of the interval for the functions 
belonging to the domain. Our task is to indicate these conditions for the doublets F 6 Dh 
at the boundaries r = and r = oo which are the so-called singular ends of the differential 
expression h ([7]), see [18]. The singularity of the left end r = is due to the nonintegrability 
of the free terms (the coefficient functions without derivatives) in h. In the case of singular 
ends, there is no universal explicit method for formulating s.a. boundary conditions. In our 
case where the coefficient function in front of the derivative d/dr is independent of r, while 
the free terms are bounded as r — > oo, the problem of s.a. boundary conditions is related 
only to the left end r = 0. 

We begin with the adjoint h* of the initial symmetric operator h^. Using the known 
distribution theory arguments (or extending the known results for scalar differential opera- 
tors [18J to the matrix differential operators), we can easily verify that the adjoint operator 
h is given by 

{ J F : F are absolutely continuous in (0, oo) , 1 
Dh * = D * = \ F,hF = Ge £ 2 (0,oc), j ' (37) 

h*F(r) = hF(r) , 

i.e., the adjoint h* is associated with the same differential expression h, but defined on a 
more wide domain Z)*, D h (p) C D*, which is the so-called natural domain for the differential 
expression h. Because the coefficient functions of the differential expression h are real, it 
follows that the deficiency indices of the initial symmetric operator h^&re equal and therefore 
s.a. extensions of do exist for any values of parameters x and q. 

It is convenient to introduce a quadratic asymmetry form A* for h by 

A, (F) = (f, h*F\ - (h*F, f\ = (f, h*F^j = 

poo poo 

= / drF + (r) (hF) (r) - dr (hF) + (r) F (r) . (38) 
Jo Jo 

The quantity A* (F) is evidently pure imaginarj@. The form A* yields a measure of the 
asymmetricity of the operator h*, it shows to what extent the operator h is nonsymmetric. 
If A* = 0, the operator h* is symmetric and therefore s.a.. This also means that is 



The quadratic asymmetry form A* is a restriction to the diagonal of the sesquilinear anti-Hermitian 



16 



essentially s.a. and its unique s.a. extension is its closure, h = which coincides with the 
adjoint, h = h* = (jij . If A* ^ 0, a s.a. operator h = (ji^j is constructed as a restriction 

of the operator h* to the domain Dh C -D* such that the restriction of A* to Dh vanishes 
and Dh is a maximum domain [22]. 

Using integrating by parts in r.h.s. in fl38l) . it is easy to verify that A* is represented as 

A, (F) = [F] (oo) - [F] (0) , (39) 

where 

[F] (oo) = lim [F] (r) and [F] (0) = lim [F] (r) 

r— »oo r— »0 

are the corresponding boundary values of the quadratic local form [F] (r) defined bj@ 

[F] (r) = -iF + (r)a 2 F(r) = -[W)g(r) -^)/(r)] = -2i9f/^)^(r). 

These boundary values certainly exist because of the existence of the integrals in r.h.s. in 
([551) : for brevity, we call them the boundary forms at infinity and at the origin respectively. 
The asymmetry form A* is thus simply determined by the boundary formcj [F] (oo) and 
[F] (0). 

We prove that for any F G -D*, we have 

lim F (r) = 0. 

r— >oo 

We first note that FgD, implies that /iF = G is square integrable together with F. It in 
turn follows that 



dF(r) ( k o o 9 



^ F(r) +io 2 G(r) 



dr \ r r 

is square integrable at infinity o It now remains to refer to the assertion that if an absolutely 
continuous F(r) is square integrable at infinity together with its derivative dF(r)/dr, then 
F(r) — > as r —>■ oo; this assertion is an evident generalization of a similar assertion for 
scalar functions. Therefore, the boundary form at infinity is identically zero: for anyF 6 D*, 
we have 

[F] (oo) = 0, 



asymmetry form lj* defined by 

u*{F u F 2 ) = (F u h;F 2 ) - {hlF u F 2 ) 



drF+ (r) (hF 2 ) (r) - / dr {hF x y (r) F 2 (r) , 
o Jo 

the forms A* and W* define each other [22] , 

19 I.e., a domain that does not allow further extending with the condition A» = conserved. 

20 The quadratic local form [-F](r) is a restriction to the diagonal of yhe sesquilinear anti-Hermitian local 
form 

[Fi,F 2 ](r) = —iFa 2 F 2 (r) = —[fi{r)g 2 (r) — gi(r)f 2 (r)]. These two forms define each other. 
21 The representation (|39| is a particular case of the so-called Lagrange identity in the integral form. 
22 A doublet F(r) is square integrable at infinity if J R drF + (r)F(r) < oo for sufficiently large R. 
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and the asymmetry form A* is determined by the boundary form at the origin: for any 
F G -D*, we have 

A* (F) = - [F] (0) = (fg-gf)\ r=0 . (40) 

For evaluating this boundary form, we have to find the asymptotic behavior of the doublets 
F G -D* at the origin. It turns out that the doublets F can tend to zero, be finite, or infinite, 
or even have no limit as r — > depending on the values of the parameters x and q: at fixed 
j, we must distinguish four regions of the charge q that are defined by the two characteristic 
values q U j and q C j (q U j < q C j) of the charge, 



T = 7 = - Z uj = 137V7(JTT), 

q cj = |x| = j + X - ^ T = <{=► Z cj = 1370' + 1/2). 

The evaluation of the asymptotic behavior of F G -D* at the origin is based on the 
following observation. According to definition (1371) . the doublets FGfl, can be considered 
as square-integrable solutions of the inhomogeneous differential equation 

hF (r) = (-ia 2 4- + -o x - - + ma 3 ] F (r) = G (r) , (41) 
\ dr r r J 

with the r.h.s. G belonging to £ 2 (0, oo) and therefore locally integrable, which allows apply- 
ing the general theory of differential equations (see e.g. [18]) to equation (|4ip . For estimating 
the asymptotic behavior of F (r) at the origin, it is convenient to represent (j4ip as 

h_F (r) = G_ (r) , (42) 

where 

= -ta 2 d r + ^a 1 - - , GL (r) = G (r) - ma 3 F (r) G £ 2 (0, oo) . 

Let Ui and be linearly independent solutions of the homogeneous differential equation 
h-U = 0, 

U\(r) = (mr) r+ u + , q > 0, 

f (mr)- T + M _, g>0, g ^ g CJ (T+ ^ 0) , 
W - j M (o» (r) = ( InW - C/ fe j _ , = fe (T+ = 0) . (43) 

Any solution F (r) of inhomogeneous differential equation ( 1421) can be represented as 

F(r) = Cl Ui{r) + c 2 U 2 {r) + h(r) + J 2 (r) , (44) 
where C\ and c 2 are some constants and 

f 2Y7/;°[f/i(r)®f/ 2 (l/)]G_(y)^, 0<g<g nj , 

J r [C/i(r) ® U 2 (y)]G-(y)dy, g > g uj , g ^ g cj , 
I 9tf Jo PMO ® U 2 {y)]G-{y)dy, q = q cj , 
t ( r \ = f 2T+ JlTt^aCO <8> Uxiy^G^dy , g > 0, g ^ g ci , 
21 J I / r [^2(r) ® U^G-tody, q = q cj , 1 ° J 
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where <8> is the symbol of tensor product, such that [U%(r) ® ^2(2/)] is a 2 x 2 matrix , and 
r > is a constant. It turns out that the boundary form [F] (0) is determined by the two 
first terms in r.h.s. in representation (1441) and essentially depends on the parameter T. 

4.2 First noncritical region 

The first noncritical region of the charge is defined by the condition 

< q < q u j T + = 7 > ^ . 

4.2.1 Self-adjoint radial Hamiltonians 

The representation given by formulas (T4"3l) - (|45p allows evaluating the asymptotic behavior 
of F G D* at the origin. According to fj43l) . the doublet U\{r) ~ r 7 is square integrable at the 
origin, whereas the doublet ^(r) ~ r -7 is not. Using the Cauchy-Bounjakowsky inequality 
for estimating the integrals Ii(r) and hif) ( |45T) . we find 

I^r) = 0(r 1/2 ) , I 2 (r) = 0(r 1/2 ), r -> 0. (46) 

It follows that for F (r) to belong to the space £ 2 (0, 00), it is necessary that the coefficient 
C2 in front of U2 (r) in (1441) be zero, c 2 = 0, which yields 

F(r) = CiU x (r) + h(r) + 7 2 (r) = 0(r 1/2 ) -> 0, r -»■ 0, (47) 

whence it follows that for any F G -D*, we have 

[F] (0) = 0. 

This means that in the first noncritical charge region, < q < q U j, or 7 > 1/2, the operator 
h* = h is s.a. and it is a unique s.a. operator associated with s.a. differential expression h 
dZjH such that 

{ ( F : F are absolutely continuous in (0, 00) , 1 
h = \ F, hFeC 2 (0,oo), J (48) 

/iF (r) = hF (r) . 

We note that this result actually justifies the standard treatment of the Dirac Hamiltonian 
with q < v3/2, or Z < 118, in the physical literature where the natural domain for h is 
implicitly assumed 2 "*!. 

We now proceed to the spectral analysis of the obtained Hamiltonian in accordance with 
the scheme described in Subsec 14.11 A necessary short information on each item of this 
scheme is given in Appendix 14.5.21 

23 We point out a particular corollary: the deficiency indices of the initial symmetric operator are 
(0,0) in the charge region < q < q u j (7 > 1/2). 

24 The uniqueness of the Hamiltonian also implies that the notion of 5 potential for a relativistic Dirac 
particle cannot be introduced, which maybe manifests the nonrenormalizability of the four-fermion interac- 
tion. 
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4.2.2 Spectral analysis 

Guiding functional In accordance with the requirements in Appendix 14. 5. 2\ for the dou- 
blet U defining guiding functional <&(F; W) ([I]), we choose the doublet U(i)(r;W) given 
by© - (HZD, 

U(r;W) = U (1) (r;W), 

the doublet U (r; W) is real-entire,see Sec. |3j 

For V, we choose the set of doublets F(r) G Dh = -D* with a compact support. It is 
evident that V is dense in £ 2 (0, oo). 

The guiding functional $ with these U and D is simple, i.e., satisfies the properties l)-3) 
presented in Appendix I4.5.2L The property 1) is evident, the property 3) is easily verified 
by integrating by parts, and it remains to verify the property 2): the equation 

(h-E )*(r) = F (r), 

where F G T> and satisfies the condition 

$(F ;£ )= / U(r;E )F (r)dr = Q, 
Jo 

has a solution belonging to T>. 

At this point, our exposition is divided into two parts because it is convenient to consider 
the cases of 7 7^ n/2 and 7 = n/2, n = 1,2,..., separately. We first consider the case of 
7 ^ n/2, after which the extension of the obtained results to the case of 7 = n/2 becomes 
evident. 

In the case of 7 7^ n/2, any solution of the inhomogeneous equation allows the represen- 
tation 

tf( r ) = Cl U(r; E ) + c 2 U (2) (r; E ) + |- J [U(r; E ) ® U (2) (y; E )}F (y)dy+ 

+ f [ r {U i2) (r;E )®U(y;E )}F (y)dy, (49) 
A l Jo 



where U^(r; W) is given by (TLB]) - (1201) . This representation is a copy of representation (143!) 
- (H5T) for a solution of Eq. (|4"21 . where we can take r = 00 because of the compactness 
of the support of F . The integral terms in this representation have a compact support: if 
suppF C [a, b], < a < b < 00, they vanish for r > b. Choosing c\ = c 2 = 0, we obtain a 
particular solution \l/ with a compact support that has the form 

= |- y°V(r; ^0) ® £/ (2 )(y; ^0)]^^)^+ 

+ 7T / [^(2)(r;^)®?7(i/;^)]Fo(j/)dy. (50) 

Taking the asymptotic behavior of the doublets U = Um and Up) at the origin (see (1T71) . 
(|T9l ). and estimate (]4"T1) for F into account, we find that the asymptotic behavior of this 
solution as r — » is of the form ^(r) = 0(r 5 ), 5 = min(7,3/2), whence it follows that 
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Green's function To find Green's function G(r,r';W), ^ 0, of the s.a. operator h 
associated with the s.a differential expression h is to represent a unique solution \l/(r) G D h 
of the differential equation 

(h-W)V(r)=F(r) (51) 
with any F(r) G £ 2 (0, oo) in the integral form 



3C 



*(r)= / G(r,r';W)F(r')dr'. (52) 



For our purposes, it is sufficient to consider the case of > 0. As any solution, \I> allows 
the representation 

1 r°° 

*( r ) = Cl f/(r; W) + c 2 F(r; W) + — — / [U(r; W) ® 

+ ^7i7T [ r V(r;W)®U(y;W)F(y)dy, 
uj{W) J 

where V(r;W) = V(i)(r;W / ) and u(W) are given by the respective formulas ( J2~li) . ( [221 . 
and ( 1231) . This representation is a copy of representation ( l4~9l) with the change of U^) to 
V(i). It is correct because V^(r; W) with > decreases exponentially as r — > oo. The 
condition \l> G £ 2 (0, oo), which is sufficient for \l> to belong to D h (because then automatically 
h^/ = Wty + F G £ 2 (0, oo)) implies that C\ = c 2 = 0: otherwise, \I/ is non-square-integrable 
at infinity (if c\ ^ 0) or at the origin (if c 2 7^ 0) because U(r; W) with > exponentially 
grows as r —>■ 00 and V(r; W) is non-square-integrable at the origin. We thus obtain that 
the solution \1/ G -D^ of Eq. (}5T|) with any F G -D^ is represented as 



*(r) = -J— / p/(r; W) ® V{y; W))F{y)dy+ 



which is the required representation (1521) with 



-V(r;W) ®C/(r';W), r > r' 



v ' ' 7 I ^?y^(r; W 7 ) ® F(r'; W), r < r' v ' 

This expression for Green's function allows evaluating the spectral function u{E) of the 
radial Hamiltonian h and writing the inversion formulas in accordance with the instructions 
in Appendix 14.5.21 (see formulas (j2J) - (Cj)). 

Spectral function and inversion formulas According to ([7]), (1531) . and (122]) we obtain 
that 

M(c; W) = ^f{ C] W) ® V(c; W) = 

1 -U(c; W) ® U(c; W) + ^-U(c; W) ® U (2) (c; W), 



u(W) v ' v ' 2 7 
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and because U(c; E) = Z7(i)(c; E) and t/(2)(c; E) are real, formulas ([SD and ([2D then yield 

— ^ = -lim3- (54) 

for the spectral function a(E) of the radial Hamiltonian h, where lim in (|54|) is understood 
in a distribution theoretic sense as well as da(E)/dE. 

The spectral function is thus determined by the (generalized) function Qw -1 (i?), 



uT A (.E) = lim' 



u(E + is) 

At the points where the function u(E), 

uj{E) = \imu(E + is), 

e— »0 

is different from zero, we have u _1 (E) = l/u(E). 

The explicit form of u(W) f[2"3j) shows that u(E) exists and is qualitatively different in 
the two energy regions \E\ > m and \E\ < m. Therefore, we naturally distinguish these two 
energy regions in the subsequent analysis. 

We first consider the region \E\ > m. 

A direct verification shows that in this energy region, u(E) is continuous, different from 
zero, and takes complex values. It follows that for \E\ > m, the spectral function u{E) is 
absolutely continuous and 

—fir = = Q ( E » E > m > 

dE 7i uj{E) 

2 7 r(2 7 )e^r(-7 + q\E\/ik)[(x + 7 )efc + iq(E - m)](2k / m)-^ 
^ } ~ qY(-2 1 )Y( 1 + q\E\/ik)\(x- 1 )ek + iq{E-m)} ' { ' 



EI\E\, k = VE 2 - 



m 



2 



We now consider the case of \E\ < m. 
In this energy region, we have 

2 7 r(2 7 )r(- 7 - qE/r)[q(m - E) - (x + 7 )r](2r/m)-^ 
1 J ?r(-27)r( 7 - g E/r)[g(m-E)-(x-7)T] ' 1 J 

r = V m 2 — i? 2 , 

o>(i?) is real, and lim e ^ [l/u;(i? + ie)} can be complex only at the points where u(E) = 0. 
Because Y(x) does not vanish for real x, u(E) can vanish only at the points satisfying one 
of the two conditions: 

i) q(m — E) — (x + 7 )r = or ii) 7 — qE/r = —n, n = 0, 1, these are the points where 
|T( 7 — qE/r)\ = oo. 

The case i) yields E = — 7 m/x for £ = 1, but at this point, we also have — 7 — qE/r = 0, 
such that the product T(— 7 — qE/r)[q(m — E) — (x + 7 )r] ^ and cu(i?) ^ 0. 

The case ii) yields E = E n = m/ \/l + q 2 / (n + 7 ) 2 , n — 0, 1, but for £ = lat the point 
E = E , we also have g(m — E) — (x — 7 )r = 0, and consequently, |T( 7 — qE/r)[q(m — 
E) — (x- 7 )r]| < oo. 
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We thus obtain that u(E) vanishes at the discrete points 

1 (n+ 7 )2 

which form the well-known discrete spectrum of bound states. We note that the discrete 
spectrum accumulates at the point E — m, and its asymptotic form as n — > oo is 

mq 2 

e n = m-E n = —, 

which is the well-known nonrelativistic formula for bound state energies. 
In the vicinity of these points, we have 

1 Ql +0(1), Q 2 n = Jim En ~ E 



uj(E + ie) E-E n + ie w ' ^ n 

It follows that for \E\ < m, the spectral function <j(E) is a jump function with the jumps 
Q\ located at the points E = E n (the discrete energy eigenvalues (15T1) ) and 

We finally obtain that the spectrum Spec/i of the operator h is the union of the discrete 
spectrum U n {E n } C (— m, m) and the continuous spectrum containing the positive part 
[m, oo) and the negative part (— oo,m), 

Spech = (-oo, -m] U (U n {E n }) U [m, oo). (59) 

We introduce a notation 

norml ' J 1 g n t/(r; K), S = En, \E\ < m, [ ™> 
The inversion formulas ([2]), ([3]) and Parseval equality (j4]) then become 

/■oo 

y>(£) = / U novm (r; E)F{r)dr, E e (-oo, -m] U (U n {£„}) U [m, oo), (62) 

/— m /*oo 
dEU norm (r; E)ip(E) + ^Wr; E n )<p(E n ) + / dEU noim (r; E)ip(E), (63) 
-oo m 
/*oo p—ra poo 

/ \F(r)\ 2 dr= \< P (E)\ 2 dE + J2ME n )\ 2 + \v(E)\ 2 dE. (64) 

t/O «/ — oo „ «/ m 



n 



1,2,..., C = l 
0,1,2,..., ( = -1 
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Inversion formulas and Parseval equality (1621) - (1641) are conventionally treated as the for- 
mulas for the generalized Fourier expansion of doublets F e £ 2 (0, oo) with respect to the 
complete orthonormalized set of the eigenfunctions U norm (r; E) of s.a. radial Hamiltonian h 
fj48l) associated with s.a. differential expression h ([7]). 

The obtained results for the energy spectrum and (generalized) eigenfunctions coincide 
with the results obtained by the standard method based on the physical arguments: the 
energy eigenstates must be locally square integrable solutions of the differential equation 
hF = EF, their moduli must be bounded at infinity, the eigenvalues E corresponding to 
the square-integrable bound eigenstates form the discrete energy spectrum, and the non- 
square-integrable eigenstates corresponding to the continuous energy spectrum must allow 
a "normalization to 5-function" . 

As the first example, we apply these considerations to the energy region \E\ < m. In 
this energy region, the solutions of the differential equation hF = EF either exponentially 
grow or exponentially decrease (in addition,they can be non-square-integrable at the origin). 
Because the required solutions must be locally square integrable, the energy eigenstates must 
belong to £ 2 (0, oo). It is convenient to first find the solutions square integrable at infinity. 
They are given by fl2TT)-(123l). 

U (1) (r;E) + ^u;(E)U (2) (r;E) 



F{r)=cV {l) {r-E) = c 



c is a constant. These functions are square integrable at the origin and therefore on the whole 
semiaxis only under the condition uj(E) = 0, which reproduces the above results concerning 
the discrete spectrum and the corresponding eigenfunctions. 

As another example of an illustration of the standard method, we show that by a direct 
calculation of the corresponding integrals, we can establish the orthonormality relations for 
the eigenfunctions that are conventionally represented in the physical literature as 



poo 

/ U noim (r; E n )U novm (r; E')dr = 0, (65) 
Jo 

POO 

/ U norm (r; E)U nmm {r; E')dr = 5{E - E'), \E\, \E'\ > m. 
Jo 

The method for calculating is presented in Appendix 14.5.21 where it is demonstrated by the 
example of the second noncritical charge region. Unfortunately, we are unable to establish 
the the completeness relation for the eigenfunctions that is conventionally written as 

/— m 
U novm (r; E) <g> U noim (r'; E)dE + ^ ^ U noim (r; E n ) <g> U noim (r' ; E n ) + 
■°° n 

POO 

+ / U noim {r;E)®U noim {r';E)dE = 6{r-r')I, 

J m 



I is the identity 2x2 matrix, by a direct calculation of the corresponding integrals, and we 
know no heuristic physical arguments in support of the validity of this relation. 
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It now remains to consider the exceptional cases of 7 = n/2, n = 1, 2, .... As follows from 
Sec. El in a neighborhood of each point 7 = n/2 and at the point itself, we can equivalently 
use the doublets U n / 2 ) and V n (y\ (12"H1) with the change of u{W) to u n (W) ( see formulas 
( 1261) and (1281) ) and obtain exactly the same conclusions about the guiding functional and 
the same results for the Green's function, spectral function and eigenfunctions as those for 
the case of 7 7^ n/2. This evidently follows from relations ( |29l) . ( 1271) . and ( l30l) . where, in 
particular, the term ^T(— 27)v4 n (W) in the r.h.s. in (15Uj) is real for real = E. It also 
follows from these formulas that Green's function and the spectral function are continuous 
in 7 at each point 7 = n/2. 

4.3 Second noncritical region 

This domain is characterized by the condition 

q uj <q< q cj < T+ = 7 < ~. 

4.3.1 Self-adjoint radial Hamiltonians 

As in the previous section, we first evaluate asymmetry form A*(.F) (T40l) evaluating the 
asymptotic behavior of the doublets F G Dh* at the origin with the use of representation 
( BHD ~~ (ESI) . In the case ofO < 7 < 1/2 under consideration, the both U±(r) ~ r 7 and 
^2(0 ~ r ~ 7 are square integrable at the origin and estimates ( [461) hold true, such that for 
any F G D h *, we have 

F(r) = Ci(mr) 7 w + + c 2 (mr)~ 1 U- + (^(r 1 / 2 ), r — > 0, 




which in turn yields 

A*(F) = — (c^ Cl -crc 2 ). (66) 

The asymmetry form A*(.F) thus turns out to be a nontrivial ant i- Hermit ian quadratic form 
in the asymptotic coefficients c\ and c 2 , which means that operator h* (l37j) is not symmetric 
and the problem of constructing nontrivial s.a. extensions of the initial symmetric operator 
h<® (ED arises. 

In solving this problem, we follow a method in [22] that comprises two steps. 

1. Reducing the quadratic anti-Hermitian form A* as a form in boundary values or 
asymptotic coefficients c a , a = 1, 2, to a canonical diagonal form by a linear transformation 
of the coefficients c a to coefficients c + k , k — 1, ...,m + , and c_;, Z = 1, ...,m_, such that A* 
becomes 

A* = ik(^2 \c +k | 2 - ^2 l c -i| 2 )' 

1 1 

where k is some real coefficient. 
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2. Relating c + k and c_/ by a unitary m x m matrix U, 

rn 

c-i = ^c +fc , Z = 1, m, 
i 

if the inertia indices m + and m_ of the fornix are equal, m + = m_ = m. Each such a 
relation with a fixed {7 convert the form A* to zero and yields s.a. (asymptotic) boundary 
conditions specifying a s.a. extension of the initial symmetric operator, different U define 
different s.a. extensions. Conversely, any s.a. extension is specified by some U, and when U 
runs over the group U(m), we obtain the whole m 2 -parameter [/(m) -family of all possible 
s.a. extensions. 

We apply this method to our case. 

By a linear transformation 

Cl,2 — > C± = Ci ± ZC 2 , 

the asymmetry form A* is reduced to a canonical diagonal form: 

A,(F) = ^(|c + | 2 -|c_| 2 ). 
1 

Its inertia indices are (1, 1), which in particular means that the deficiency indices of with 
< 7 < 1/2 are (1,1). 
The relation 

c_ = e i9 c + , < 9 < 2tt, ~ 2tt, (67) 

with any fixed # yields boundary conditions specifying a s.a. extension hg of the operator 
h^°\ Different #'s are assigned different s.a. extensions, except equivalent cases of 6 = and 
6 = 2tt. When 9 runs over a circle, we obtain the whole one-parameter {7(l)-family of all 
s.a. extensions of the operator h^°\ 

Relation (1671) is equivalent to the relation 

Q 

c 2 = £ci, — oo < ^ = — tan - < +oo, — oo ~ +oo, 

the values ^ = ±oo are equivalent and mean that c\ = 0; we will say that ^ = oo in these 
cases. 

We let fi£ redenote the corresponding s.a. operator, = hg, and let denote its 
domain. The final result in a more extended form is formulated as follows. In the second 
noncritical region < 7 < 1/2, we have a one-parameter [/(l)-family {h^} of s.a. operators 
associated with s.a. differential expression h (jTj). They are specified by s.a. boundary 
conditions and are given by 



hp 



F(r) : F(r) is absolutely continuous in (0, 00), F, hF C C 2 (0, 00), 
Dt- = { F(r) = c[(mr) 7 w + + £(mr)~ 7 w_] + Oir 1 ' 2 ), r -»• 0, -00 < ^ < +00, 

F(r) = c(mr)^n_ + Oir 1 / 2 ), r -> 0, ^ = 00, 
ht.F = hF, 



(66 
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The inertia indices coincide with the deficiency indices of an initial symmetric operator. 
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where c is an arbitrary complex number. 

In other words, the only s.a. differential expression h does not uniquely define a s.a. 
operator in the charge region < 7 < 1/2, and an additional specification of the domain in 
terms of s.a asymptotic boundary conditions involving one real parameter £ is required. 

4.3.2 Spectral analysis 

The spectral analysis in this charge region is quite similar to the analysis performed in 
Subsec. 14.21 related to the first noncritical region 26 !. We therefore only point out necessary 
modifications and formulate the final results. 

In the case of £ = 0, the the corresponding analysis is identical to that in the previous 
Subsec. 14.21 and the results obtained there are directly extended to the region < 7 < 1/2 
and are given by the same formulas. 

Until said otherwise, we assume that < |£| < 00, £ is arbitrary , but fixed. The case of 
£ = 00 is considered separately below. 

For the doublet U(r; W) defining guiding functional dTJ , we choose the doublet 

U^r; W) = U {1) (r; W) + £[/ (2) (r; W) (69) 

satisfying the condition 

£/f(r; W) = (mr) 1 u + + ((mr) _7 u_ + 0(r~ 7+1 ), r — > 0, 

where and Up) are given by formulas ([TBI) - (1201) . As before, U^(r; W) is real-entire in 
W . The corresponding guiding functional is denoted by $z(F, W). 

For V, we choose the set T>^ of doublets belonging to and having a compact support. 

The guiding functional with the chosen U% and T>^ is simple. Indeed, the properties 
1) and 3) are evident; as to the property 2), the solution \I> of the inhomogeneous equation 
(h — Eo)^ = Fq, F G T>£, with the property \I/ G T>^ is given by a copy of (J50|) . where 
the solutions U = Un\ and Up) of the homogeneous equation are replaced by the respective 
solutions U% and Um with the Wronskian Wr(L^, Un\) = 

*(0 = J°°[U^r; E ) U {1) (y; E )}F (y)dy+ 

"A / U {1) (r;E )®Ut(y;E )]F (y)dy. 
^7s Jo 

Green's function G^(r, r'\ W) , ^sW > 0, of the operator is defined as the kernel of the 
integral representation for any \l/ G in terms of the doublet F = (h^ — W)ty G £ 2 (0, 00): 



POO 

*(r) = / G^r,r';W)F(r')dr'. 
Jo 



This representation is a copy of formula fl52|) . The natural difference is the change of U — Urn 
to U — U^ because the condition ^ G D% implies that \I> satisfies s.a. asymptotic boundary 



3 A simplifying thing is that the particular cases of 7 = and 7 = 1/2 are excluded. 
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conditions (1681) . The final result is 

Gz(r,r';W) = 



-^V(r;W)®U ( (r';W), r>r>, 



y ) U^r-W)®V{r'-W), r<r', 



where 



V(r; W) = %(r; W) = U^r; W) + ±^{W)U {2) {r; W), 

27f 

uz(W) = - Wr(U, V) = uo(W) -, 



(70) 



V(i)(r; W) and u(W) are given by ( 12TI) - ( |23|) . which is a copy of representation ( 1531) with 
the natural change of U = C/(i) and a; to the respective U = and co^. 
It follows, see P - dZJ), that 

M(c; W) = ^y^U^c; W) ® V(c; W) = 

' : £^(c; W) ® t/ £ (c; W) + ^-U {2) {c- W) ® C/ C (c; W), 



?v ' w ?v ' ' 2 7 

and, because the both L^(c; i?) and t/( 2 )(c; i?) are real, that the spectral function u^(E) of 
the radial Hamiltonian h^, < |£| < oo, is given by 

dE 7r £-*o LO^E + ie)' 
a copy of expression (1541) with the change of to cj^. The spectral function is determined 



by the (generalized) function Sco , £ 1 (£'), 



cj £ 1 (E) = lim — — — = lim 



? e->o ll>s(E + ie) e^ouj(E + ie) 

At the points where the function 

ujz(E) = \imu ( (E + ie) = \im[u(E + ie) - 2^/q] 

is different from zero, we have (E) = l/cu^E). Because to^(E) differs from u(E) by a real 
constant — 2^^/q, the two energy regions > mand \E\ < mare naturally distinguished as 
before, and the corresponding analysis in each region is similar to the analysis in the previous 
subsection. In the energy region \E\ > m, the function uj^(E) is continuous, different from 
zero, and complex as well as u(E) ( 1551) . The spectral function cr^{E) for \E\ > m is therefore 
absolutely continuous, and 

d ^ {E) - If* 1 -C*(E) 

which is an analogue of fl55|) with the change of ui{E) and Q(E) to the respective uj^(E) and 
Qt(E). 
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In the energy region \E\ < m, function u(E) (I56p is real, and therefore, the function 
uj^(E) is also real. As in the case of the first noncritical charge region, it follows that for 
\E\ < m, the spectral function cr^(E) is a jump function with the jumps Q 2 n located at the 
points E = E^ n) the discrete energy eigenvalues, where ui^{E^ n ) = 0, and 



Q £n = hm 



UJ^E) 



As a result, we obtain that 

da^E) 
dE 



J2QlJ(E-E^ n ), \E\< 



m, 



which is an analogue of ( 1581) with the change of E n and Q n to the respective E^ n and Q^ n - 
Unfortunately, we are unable to find an explicit formula for the discrete energy eigenvalues 
E^ n with £ ^ 0, we only note that, as in the first noncritical charge region, there are infinitely 
many of such levels accumulating at the point E = m, and their asymptotic form as n — > oo 
is given by the previous nonrelativistic expression independent of £ : 

_ mq 2 
£{,n = m- E^ n = —. 

The lower bound state energy essentially depends on £, and there exists a value of £ 
for which the lower bound state energy coincides with the boundary E = —m of the lower 
(positron) continuous spectrum. 

The whole spectrum Spec/i^ of the radial Hamiltonian is given by a copy of (|59p with 
the change of E n to E^ n . 

The inversion formulas and the Parseval equality 

POO 

tpt(E) = / t/ ? ,norm(r; E)F(r)dr, E e (-oo, -m] U (\J n {E tn }) U [m, oo), 
Jo 

/— m 
dEU 6>noim (r; E)(p^{E) + ^ U 6>n0Im (r; E^ n )^(E^ n ) + 
■°° n 

poo 

,iionn 

(r; E)(p^(E), 

J m 

poo p—m poo 

/ \F(r)\ a dr= \^(E)\ 2 dE + J2\^(Es,n)\ 2 + \^(E)\ 2 dE, 

Jo J — oo J m 

are written in terms of the normalized eigenfunctions U^ noTm (r; E) and the Fourier coefficients 
f^(E) that are defined by copies of formulas (I60I) and (l6Tj) with the addition of the subscript 
£ to all the symbols Q, U, Q n , E n , and $. These relations are copies of formulas (!62|) -( l64|) . 

As before, the energy spectrum and (generalized) eigenfunctions of the radial Hamiltonian 
fi£ can be obtained by the standard method using the physical arguments. 

As an example, we consider the energy region \E\ < m where the solutions F of the 
differential equation hF = EF either exponentially grow or exponentially decrease as r — ► oo, 
all solutions being square integrable at the origin. Only exponentially decreasing solutions 
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F(r) = cV (1) (r; £) = c[C/ (1) (r; £) + |-o;(£;)C/ (2) (r; £)] = 
= c[^(r;E) + |-^( J E)f/ (2) (r; J E)], 

c is a constant, are proper. It is remarkable that such solutions are square-integrable on 
the whole semiaxis for any energy values E G (— m, m). But they satisfy s.a. asymptotic 
boundary conditions 0681) only if uj^(E) = 0, which reproduces the results for the eigenval- 
ues and eigenfunctions of the discrete spectrum of the operator h^. We note that the s.a. 
boundary conditions have a physical meaning of the condition that probability flux density 
vanish at the boundary, the origin in our case . We did not refer to this requirement in the 
first noncritical charge region because it is automatically satisfied in this region. 

We can also establish the orthonormality relations for the eigenfunctions U^ noTm (r;E), 
which are copies of relations fl65l) . by a direct calculation of the corresponding integrals. This 
calculation illustrating the general method applicable to all the charge regions is presented 
in Appendix 14.5.21 

As before, we are unable to establish the completeness relation for the eigenfunctions by 
a direct calculation or based on heuristic physical arguments. 

We now touch briefly on the case of £ = oo where the s.a. asymptotic boundary conditions 
for any F G .D^are 

F(r) = c(mr)~ 7 u_ + 0(r 1/2 ), r -> 0. 
For the doublets U and V, we choose the respective 

U(r;W) = U 00 (r;W) = U {2) (r;W) } 
U^ir- W) = (mr)" 7 u_ + O^ 1 " 7 ), r -> 0, 



and 



V(r; W) = J^-VfoCr; W) = U {2) (r; W) - l-u^U^r, W), 
u) 0B {W) = -Wr{U,V)= l ~ r ' 



q 2 u(W)' 



Performing calculations completely similarly to those in the case of |£| < oo, we find the 
spectral function a^E), 

— = — limb — r = hm isulE + is). 

dE n e^o Uoc(E + ie) 7r47 2 e-o 

All other results concerning the structure of the spectrum and inversion formulas are 
also completely similar to the results in the case of |£| < oo . In particular, the bound state 
spectrum is determined by the poles of the function uj(E) in the energy region \E\ < m; it 
can be evaluated explicitly. 

In conclusion, we note that the spectrum and the normalized eigenfunctions U^ norm (r; E) 
are continuous in £ , the points £ = and £ = oo included. 
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4.4 Critical charges 

This region is defined by the charge values 

q = Qcj = | x | — j + 1/2 T + = 7 = 0. 

The charge values q = q C j stand out because for q > q C j, the standard formula (I5"7j) for the 
bound state spectrum ceased to be true yielding complex energy values. But we will see 
that from the mathematical standpoint, no extraordinary happens with the system for the 
charge values q > q C j , at least, in comparison with the previous case of q U j < q < q C j. 

4.4.1 Self-adjoint radial Hamiltonians 

In constructing the s.a. Hamiltonians with critical charges, we literally follow the method 
presented in the previous Subsec. 14.31 where the second noncritical charge region was consid- 
ered and therefore restrict ourselves to only key remarks. 

It follows from representation fj43l - (145]) that the asymptotic behavior of the doublets 
F G D* in the case of q = q C j is given by 

F(r) = Cl u + + c 2 u^\r) + 0(r 1/2 In r), r -> 0, VF G D h *, 
which yields the expression 

A*(F) = — (crc 2 -cid) 

for A*(F), this expression is completely similar to the previous case, see fl66l) with the changes 
27/g -> l/q cj and c x c 2 . 

Therefore, in the case of critical charges, i.e., for 7 = 0, we also have the one-parameter 
[/(l)-family {h^}, —00 < £ < +00, of s.a. operators associated with s.a. differential ex- 
pression h ([7]) and specified by s.a. asymptotic boundary conditions, 




(F(r) : F(r) is absolutely continuous in (0, 00), F, hF C £ 2 (0, 00), 
F(r) = c{u m {r) + £u + ) +0(r 1 / 2 lnr), r -»• 0, -00 < ^ < +00, 
F(r) = cm + + 0{r x l 2 lnr), r — > 0, £ = 00, 



(71) 

Z)^ denotes the domain of the operator /i(o)£, and £ = 00 corresponds to the equivalent 
cases of £ = +00 and £ = —00. 

4.4.2 Spectral analysis 

The spectral analysis follows the standard way presented in the previous subsections, and 
therefore, we only cite the final results. 

We first consider the case of £ 7^ 00. For the doublet U(r; W) defining the guiding 
functional ([T]), we choose the doublet 

Uf\r; W) = U$(r; W) + £C/ (1) (r; W) 
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with the asymptotic behavior 

Uf\r]W) = u {0) (r) + £w+ + O(rlnr), r -> 0, 

where the doublets C/(i) and are given by formulas (13TT) - fl34l) ; uj?'(r; W) is real-entire 
in W. 

For T>, we choose the set T>^ of doublets belonging to and having a compact support. 

The guiding functional with these and is simple. In particular, the solution 

of the inhomogeneous equation (h — E )ty = F , F Q G X>t , is given by a copy of (J50j) 

with the change of U^,U^), and j- = — (Wr(Z7(i), f/(2))~ 1 to the respective L^ ' 1 , Z7(i), and 

-^.^-(WrCC/f,^))- 1 . 

Green's function Gf\r,r';W), > 0, of the Hamiltonian h(o)£ is defined by a copy 
of fl53|) with the change of U = Um, V = Vm, and = — Wr(/7(i), V(i)) to the respective 

17 = flf , V = V$ , and 4 0) = -Wr(C/f , ^ 
Gf (r,r';W) = / ^ 



) 

W^V^(r; W) ® Uf\r'; W), r > r', 



i rr(°)/..w^v(°) 



r;W0 <8> V ( \y(r'; W), r < r', 



see formulas (1351) and (|36|) ; is conveniently represented as 

Vg^r; = L/f >(rj W) + q CJ uf\w)U (1) (r; W), (72) 

where 

u f\W) = -Wr(u(°\V$) = J°\W) - J-, 
is given in (|36|). 

The spectral function cr^^-E 1 ) of the radial Hamiltonian h^is defined by 

d°f\E) l_ 1 



— lim 



dE vre^o ..(0) 



and is determined by the (generalized) function 1 (E), 
uf^ ~ 1 (-E) = lim —-- = lim 



^ cuf (E + is) ^° ^ (0) (E + is) - i/q cj ' 
At the points where the function 

uf{E) = \imuf{E + is) = uj^(E) - Hq c , 

with 

J°\E) = \imu {0) (E + ie), 
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is different from zero, we 



The two energy regions \E\ > mand \E\ < to are naturally distinguished as before. 
In the region \E\ > to, the function co^(E) is given by 



a;® 



(E) = -1 {ln[2e— ^k/m] + H-iq CJ \E\/k) + - 2^(1) + C/(2<fe) ) , 

Qcj l ^q C j& ) 



e = E/\E\, k = \[W 



to . 



It is continuous, different from zero, and complex, therefore, the spectral function a^(E) 
for \E\ > m is absolutely continuous, and 

daf\E) 11 , m 

In the region < to, the function c<j^(i?) is given by 

J°\E) = -1 (ln(2r/TO) + ^(-^/t) - - 2^(1) + C/(2^) ) , 

q C j I ^q C j& J 

r = V to 2 — E, 

it is real, and therefore, the function uj^\e) is also real. As in the previous cases, the 
spectral function a^(E) for \E\ < to is a jump function with the jumps 

L^, n J - urn (0) 
located at the discrete energy eigenvalues E^l where uf\Ei l) = 0, such that 



dE 



As in the previous case, we are unable to find an explicit formula for E^ n (except the 
case of £ = oo, see below), we only note that there exists an infinite number of of the discrete 
levels which accumulate at the point E = to, and their asymptotic behavior as n — > oo is 
described by the same nonrelativistic formula: 

(o) _ p (0) _ mq 2 

The lower bound state energy essentially depends on £, and there exists a value of £ 
for which the lower bound state energy coincides with the boundary E = —to of the lower 
(positron) continuous spectrum. 

All other results concerning the inversion formulas and Parseval equality are written in 
terms of the normalized (generalized) eigenfunctions U^ orm (r; E) and the Fourier coefficients 
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if ® (E) as copies of relations (j5Dj) - ([51j) with the addition of the subscript £ and superscript 
(0) to all the symbols Q, U, Q n , E n , and $. 

As before, the energy spectrum and eigenf unctions of the radial Hamiltonian hm)^ can 
be obtained by the standard method using the physical arguments. As an example, we 
again consider the energy region \E\ <m where the solutions F of the differential equation 
hF = EF either exponentially grow or exponentially decrease and are square integrable at 
the origin. The exponentially decreasing solutions are described by the doublets 

F(r) = cV$(r; E) = c[uf\r; E) + q cj cuf (E)U {1) (r; E)}, 

where c is a constant. They are square integrable, F G £ 2 (0, oo), for any E in the interval 
E G (— m, m), but satisfy s.a. asymptotic boundary conditions (17TT) only if uj^(E) = 0, 
which reproduces the results for the eigenvalues and eigenfunctions of the discrete spectrum. 
We can also verify the orthonormality relations for the eigenfunctions f^°n 0rm (^; E) that are 
analogues of relations fl65l) by a direct calculation of the corresponding integrals with the 
method described in Appendix 14.5.21 

We touch briefly on the case of £ = oo where the s.a asymptotic boundary conditions for 
any WF G are 

F(r) = cu + + 0(r 1/2 In r), r -> 0. 
For the doublets U and V, we choose the respective 

U(r; W) = u£>(r;W) = U {1) (r;W), 
U£\r; W) = cu + + O(r), r -» 0, 



and 



V(r; W) = — w ^ V$(r; W) = U$(r; W) - q^dg (W)U (1) (r; W), 
to£\W) = -Wr(U {lh V)= ' 



q 2 cj cuW(Wy 

Proceeding completely similarly to the case of |£| < oo, we find the spectral function ag(E): 

= Ii im!3 ^ = hm^^ + ze). 

dE ne-,o JVfi + ie) * 3 ^° 

The structure of the spectrum, the inversion formulas, and the orthonormality relations 
are also completely similar to the corresponding results in the case of |£| < oo . In particular, 
the bound state spectrum is determined by the poles of the function (E) in the interval 
\E\ < m. It can be evaluated explicitly and is given by formula (15TI) with 7 = 0, the energy 
of the lower level with ( = — 1 is equal to zero, E^ = 0. 

In conclusion, we note that the spectrum and the normalized eigenfunctions U^°l orm (r; E) 
are continuous in £, the point £ = oo included. 
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4.5 Overcritical charges 

This region is denned by the charge values 



Q > Qcj — M — j + 1/2 -<=>- T + = ia, a = \J q 2 — x 2 > 0. 

In constructing s.a. Hamiltonians and analyzing their spectral properties in this charge 
region, we canonically follow the methods used in the previous cases and therefore only cite 
the main results. 

4.5.1 Self-adjoint radial Hamiltonians 



According to representation (1431) - ( 1451) . the asymptotic behavior of the doublets F G -D* in 
the case of q > q C j is given by 

F(r) = c l (mr) ia u + + c 2 (mry ia U- + 0(r 1/2 ), r -> 0, VF G D h *, 
-( 1 ^1 



which yields 



A*(^) = — (|ci| 2 -|c 2 | 2 ) 



It follows that in the case of overcritical charges, i.e., for T + = ia, we have the one- 
parameter {7(l)-family {he}, < 9 < n, ~ n, of s.a. operators associated with s.a. 
differential expression h ([7]) and specified by s.a. asymptotic boundary conditions 



F(r) : F(r) is absolutely continuous in (0, oo), F, hF C £ 2 (0, oo), 
D e = { F(r) = c[e id (mr) ia u + + e~ ie {mr)- i<T u^ + 0(r 1 / 2 ), r -»• 0, 

< < 7T, ~ 7T, 



(73) 



Dg is the domain of hg. 
4.5.2 Spectral analysis 

For the doublet U(r; W) defining guiding functional (CD), we choose the doublet 

U e (r; W) = e* e U {1) (r; W) + e~ ie U m {r; W) 

with the asymptotic behavior 

U(r; W) = e id {mr) la u + + e~ ie {mr)' la u^ + 0{r), r -> 0, 

where f/(i) and f/( 2 ) are given by formulas ([TB]) - (J2"0j) (with T + = ia ); f/e(r; W 7 ) is real-entire 
in because £7(2) = £7(i) for T + = icr and real W = E. 



27 The relation C2 = e l9 ci, < 8 < 2ir, denning s.a. boundary conditions (compare with (IST)) ). is equivalent 
to the relations c\ = e c, C2 = e~ %t c, < 8 < ir, with the change 8 — > 27r — 26*. 
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For T>, we choose the set Vg of the doublets belonging to D e and having a compact 
support. 

Using the doublets Ug(r; E ) and Um(r; E ) for constructing the solution \1/ e Vg of the 
inhomogeneous equation (/i — E Q )ty = F , F G we verify that the guiding functional $0 
is simple. 

Green's function Gg(r,r'; W), > 0, of the Hamiltonian hg is constructed in terms of 
the doublets U = Ug and V = Vg, where 

Ve(r; W) = e - ie + lg~ {w) V(i)(r; W) = U e (r; W) - £ug(W)U (r- W), 
Mr; W) = -\e ie U {1) {r; W) - e^f/ (2) (r; W)\, 

2zcr g l+u>(W)e 2iy 

V(i)(r; W 7 ) and w(W) are given in (j2"T]) - ([221) (with T + = zV); ^(r; W 7 ) is real-entire in VF. 
As a result, we obtain that 

, f ^(r;W)®^(r';W),r>r', 

" 1 ^C/ 9 (r;Ty)®^(r';W), r < r'. 

The spectral function ag(E) of the radial Hamiltonian hg is defined by 



and is determined by the (generalized) function Qw -1 ^), 

wr 1 ^) = lim — — ^ -. 

£->o U g{E + is) 

At the points where the function 

LOgiE) = limug(E + is) 

e— >0 

is different from zero, we have wj\E) = l/u g (E). 

The two energy regions \E\ > mand \E\ < mare naturally distinguished as before. 

In the region \E\ > m, the function Ug{E) is continuous, different from zero, and complex, 
therefore, the spectral function <Tg(E) for \E\ > m is absolutely continuous, and 

^=I 9 -1 sfl }(B). 

dE 71 UJg(E) 

In the region \E\ < m, we have 

Q{E) = n2ia) n-ia - Eq/r) r(x + io) - g(m - E) {2r/m) -2u, s e - 2 *e(*) 
V ; r(-2z(j) T(ia-Eq/r) t(k - ia) - q(m - E) K 1 ! 
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therefore the function 

ug(E) = — tan(e(E)-9) 
Q 

is real. 

It follows that the spectral function <Jq{E) for \E\ < m is a jump function with the 
jumps Ql 



,y> y E B , n -E 



located at the discrete points Ee, n where ujg(Eg in ) = 0, such that 

^ = £«»,»^ -£«.»>- 

n 

We failed to find an explicit formula for the discrete energy eigenvalues Eg tn . We only note 
that there is an infinite number of the discrete levels which accumulate at the point E = m. 
Their asymptotic form as n — > oo is given by the previous nonrelativistic formula: 

_ mq 2 

£d,n = 111 — Eg tn = — — ^ ■. 

The lower bound state energy essentially depends on 9, and there exists a value of 9 for which 
the lower bound state energy coincides with the boundary E = —m of the lower (positron) 
continuous spectrum. 

The inversion formulas and Parseval equality in terms of normalized (generalized) eigen- 
functions Ug tIl0Tm (r; E) and the Fourier coefficients (pg(E) are copies of formulas fl60|) . (lB"Tj) . 
(JS25, fl, and ([MD. 

As in the previous subsections, a comment on the applicability of the standard method for 
finding the energy spectrum and eigenfunctions based on the physical arguments holds. As 
an example, we consider the energy region \E\ < m where the solutions F of the differential 
equation hF = EF either exponentially grow or exponentially decrease as r — > oo, any 
solution being square integrable at the origin. Only exponentially decreasing solutions 

F = cV e (r; W) = c[U e (r; W) - ±aj e (W)U e (r; W)], 

c is a constant, are proper. They are square integrable, F G £ 2 (0, oo), for any E in the 
interval \E\ < m, but satisfy s.a. asymptotic boundary conditions ( 1731) only if u>$(E) = 0, 
which reproduces the results for the eigenvalues and eigenfunctions of the discrete spectrum. 
We can also establish the orthonormality relations for the eigenfunctions t^6»,norm( r ; E), which 
are analogues of relations (165]) . by a direct calculation of the corresponding integrals with 
the method described in Appendix I4.5.2L 

In conclusion, we point out that the number of s.a. extensions of the total Dirac Hamil- 
tonian (i.e., the number of independent parameters of s.a. extensions) depends on the charge 
values q (on the value of Z) as follows: it is easy to verify that in the interval q n < q < q n +i, 
where 

"0, n = 0, 
V 7 ^ 2 - 1/4, n = l,2, 



'1 • • • 5 



37 



the number of independent parameters of s.a. extensions is equal to 2n. This follows from 
the fact that the total Dirac Hamiltonian is a direct sum of its parts unitary equivalent to 
the radial Hamiltonians, see (P), flHD- 

Appendix A 

One of the methods for finding the spectrum of a s.a. differential operator and its complete 
system of (generalized) eigenfunctions and constructing the corresponding Fourier expansion 
with respect to these eigenfunctions (the so-called inversion formulas) is based on the Krein 
method of guiding functionals. For ordinary scalar differential operators, this method is 
described in [18], but it is directly extended to the case of ordinary matrix operators. We 
here present the key points of the method as applied to our case where it suffices to consider 
only one guiding functional. This implies that the operator spectrum is simple^!; we call 
such a functional a simple guiding functional. 

By definition, a guiding functional <&(F; W) for a s. a. operator h associated with a 
differential expression h is a functional of the form 

POO 

$(F;W) = U{r;W)F(r)dr, (1) 
Jo 

where U(r\ W) is a solution of the homogeneous equation 

(h-W)U = 

which is real-entire in W and F(r) belongs to some subspace T> C D PI D^, where D is a 
space of doublets with a compact support such that V is dense in £ 2 (0, oo). 

A guiding functional $>(F; W) is called simple if it satisfies the following conditions: 

1) for a fixed F, the functional <&(F; W) is an entire function of W; 

2) if 

$(F ; E ) = 0, %E = 0, F e V, 
then the inhomogeneous equation 

(h - E )V = F 

has a solution $6D; 

3) 3>(hF;W) = W$(F;W). 

We note, that the existence of a simple guiding functional is conditioned by the existence 
of appropriate U and T>. A guiding principle for the choice of U(r; W) is that its behavior 
as r — ► must conform to the asymptotic behavior admissible for the doublets belonging to 
D h , or, roughly speaking, U(r;W) at the origin must belong to D h . This corresponds to 
the conventional physical requirement that the (generalized) eigenfunctions of the operator 
h, being generally non-square-integrable, but "normalizable to 5-function", satisfy the s.a. 
boundary conditions specifying h. 

If a simple guiding functional exists, then the s. a. operator h has the following spectral 
properties: 



For a definition, see [18j . In the physical terminology, this means that the spectrum is nondegenerate. 
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1) the spectrum of h is simple, and there exists a spectral function u(E), a nondecreasing 
real function continuous from the right and such that the set of spectrum points coincides 
with the set of growth points^) of the function cr(E); 

2) the inversion formulas 



and the Parseval equality 



$(E) = / U{r;E)F(r)dr, (2) 
Jo 

/oo 
U(r;E)$(E)da(E), (3) 
-OO 

roo /»oo 

/ \F(r)\ 2 dr= / \<&(E)\ 2 do(E) (4) 

J0 J -oo 



hold, where G L%(— oo, oo), F(r) G £ 2 (0, oo), the function U(r;E) in the integrands 

in d2J) and (J3]) can be defined as zero outside of the spectrum points of the operator h 
(outside of the growth points of a(E)), and the convergence of the integrals in <^ and (J3]) 
in general must be understood in the sense of the convergence with respect to the metrics of 
the respective spaces L 2 (— oo, oo) and £ 2 (0, oo). This means that the set of the (generalized) 
eigenfunctions {U(r; E),E G Spec/i} of the operator h forms a complete orthogonal system. 

The spectral function u{E) can be expressed via the resolvent of the operator h. As is 
known, see |18j . the resolvent R(W) = (h — W)' 1 with ^sW ^ is an integral operator with 
the kernel G(r,r'; W) ( Green's function). The spectral function a(E) is expressed in terms 
of Green's function as follows: 

U(c;E)®U(c;E)da(E)=dM(c;E), (5) 

I rE+S 

M(c;E)= lim lim - / 3M(c; E' + ie)dE', (6) 

M(c;W) = G(c-0,c + 0;W), (7) 

where c is an arbitrary internal point of the interval (0, oo). We note that for any E, one 
of the diagonal elements of the matrix Z7(c; E) <S> U(c; E) is different from zero. Of course, 
cr{E) is independent of c. 



Appendix B 

In this Appendix, we describe a method for calculating the so-called overlap integrals for 
the solutions of the differential equation hF = EF . It is based on the integral Lagrange 
identity and on evaluating the asymptotic behavior of the solutions at the boundaries, the 
origin and infinity. We apply this method to proving the orthonormality relations for the 
(generalized) eigenfunctions of the radial Hamiltonian and illustrate it by the example of the 
second noncritical charge region. 

29 The set of growth points of the function cr{E) is the complement of the open set of constancy points of 
the function <r(E). A point Eq is the constancy point of the function cr(E) if there exists a vicinity of the 
point Eq where <r(E) is constant. 
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We call the integral 



poo pit 

/ F(r; W)F'(r; W')dr = lim / F(r;W)F'(r;W')dr 

Jo i?.^oo,e-->0 J t 

for two doublets F and F ' the overlap integral for these doublets. Let F and F ' satisfy the 
respective homogeneous equations 

(h - W)F(r; W) = and (h - W')F'(r; W) = 0. 

Then the equality for the overlap integral 

poo 

/ F{r; W)F'{r; W')dr = I°° - J , (1) 
Jo 

where 

r- = u m ^Ln, (2) 

r ^oo W - W 

Wr(r; F , F') 

h = H W-W> (3) 

holds. The equality is a special case of the integral Lagrange identity. In what follows, 
we are interested in the case of real W and W, W = E and W = E' . The overlap integral 
is understood as a generalized function of E and E' . Evaluating the overlap integrals is thus 
reduced to evaluating the asymptotics of the Wronskian of the corresponding doublets at 
the boundaries. 

We begin with evaluating the asymptotics of some basic functions and doublets. 
Let \E\ > m, in which case we have K = \E\k/E, k = \/E 2 — m 2 , and let r — > oo. Using 
the known asymptotics of the functions P',x) (see, for example, [21]), we have 

$(a,/3;2^r) - — — !-^t^^ 

a = T + % 2 qE/k, & = ±1, £ 2 = ±1, = 1 + 2T, 

where T is any real or pure imaginary number, T ^ —n/2, n = 1,2,.... If T is real, 
T = 7 7^ —n/2, n = 1,2, then 

(mr) T $+(r, T, E, k) -> 2A(T, £) cos^(r; T, £), 

(mr) T $_(r, T, £, k) -> ^A(T, £) sin^fr; T, £), 



where 



r(l + 2T)(2fc/m)- T e-^/( 2fc ) 
A l 1 > - — 



\T(l + T + iqE/k)\ 

QE , . , . 7T 

h V In 2A;r - - 
k 

ip r (T, E) = argT(l + T + iqE/k) 



V»(r; T, E) = kr + ln(2Ar) - ^ - ^ r (T, £), 
A; 2 



In this Appendix, a notation like F 1 denotes another function, but not a derivative. 
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which yields 

X(r,T,E,k)^A(T,E) 



sin ?/;(>• T E) / ( m +^)(^+7) 
cos^(r;T,EK + Sm ^ r ; 1 '^ ( * 

/e \ m — E 



Let |E| < m. For our purposes, it suffices to know that the doublets U (r; E n ) and V(r; E) 
decrease exponentially as r — > oo. 
We use the relation 



£ sm[-0 r (r; T, £) ± ^ r ( r ; T ', \ E \ sm[(ifc ~ k ') r 



\E\ 



E-E' 



k 



k-k' 



^-n5(k - k') = 7r5(E - E'), r -> oo, 



(4) 



which holds in a distribution-theoretic sense. We call the expressions of the form 

a±(E, E') sii# r (r; T, E) ± ^ r (r; T', 
6±(S, £') cos[^ r (r; T, £) ± ^ r (r; T', E% 

where a±(E,E'), b±(E,E') are finite at E — E', the quickly oscillating expressions (QO); 
such expressions have the zero limit in a distribution-theoretic sense as r — > oo. 
This allows obtaining the limits J 00 for the basic doublets. We have 

U (1) = U (1) (r;i,E), Ufa = Ufa(r; 7, E'), 7 > 0, 
{E-E'Y'Wir-U^ Ufa) -> 



-> A('y,E)8(E-E'), r^oo 

U {1) = U w (r;>y,E), Ufa = Ufair^,^), < 7 < 1/2, 
(E — £'') -1 W(r; U^),Ufa) — »■ 

£ su# r ( 7 ,£)-^ r (- 7 ,£')] 



-S( 7>J B) 



|£| n(E-E') 

7/c cos[^ r (7, -E) - ^ r (-7, £")] 



7r(E - £?') 



+ QO 



-8(7, E)5(E — E') — -8(7, 



7/c cos[^ r (7, - ^ r (-7, £')] 



, T > OO, 



(5) 



(6) 



U m = U (2) (r ;i ,E), Ufa = Ufa(r; % E f ), < 7 < 1/2, 



'(2)1 
n/\-l 



(1) 



(£-£')-W(nf/( 2 ),4) 



B(^,E)5(E' - E) - B(-f,E) 
B( 1 ,E)5(E-E') + B( 1 ,E) 



7k cos[^ r (7, E') - ^ r (-7, £)] 



7/c cos[^ r (7, -E) - Vr(-7, 



r — > 00, 



(7) 
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U m = U (2) (r- n ,E), U' {2) = f/ ( ' 2) (r; 7 ,£'), < 7 < 1/2, 
(E-E')- l W{r-U {2) ,U{ 2) )^ 

^A{- ll E)5{E~E') 1 r^oc, (8) 



where 



A T, £?) = A (T, E) — yHc] s ^| 1 1 s S 9 

E(7,£) = ^A( 7 ,£)A(-7,£). (10) 

For the corresponding limits Jo, we respectively find (the relations are valid for any E, E') 

U<i) = f/(i)(r;7,E), Ufa = Ufa^^E'), 7 > 0, 

(E - E'y l W{r- U {1) , Ufa) -+ 0, r - 0, (11) 



C/ (1) = U {1) (r;j,E), U{ 2) = Ufa(r; 7, E'), < 7 < 1/2, 

(J5 - E'y'Wir; U {1) , Ufa) -> - . J* 7 - , r - 0, (12) 



17(2) = [/(2)(r; 7 ,£), l/fo = Ufa(r;j,E'), < 7 < 1/2, 

(J5 - JST 1 ^; tfa, tf a) ) - r - 0, (13) 



Up) = U {2) (r;j,E), U{ 2) = Ufa(r; 7, E'), < 7 < 1/2, 

(E - ET'Wir; U {2) , U{ 2) ) - 0, r - 0. (14) 

The obtained relations allow calculating the overlap integrals and proving the orthonor- 
mality relations for the eigenfunctions of the radial Hamiltonians. As an example, we con- 
sider the second noncritical charge region, the other charge regions, including the critical 
and overcritical regions, are considered quite similarly. 

We have to calculate the integral 



Uz(r;E)Ut:(r;E')dr, 

where U^(r; E) is defined by Eq. ( |69|) . Using relations (JT|) (J3j) and §E§ - ( |T4l) . we find 

POO 

/ C/ C (r; E)Uz(t\ E')dr = C^(E)5(E - E'), \E\, \E'\ > m, 
Jo 

C^E) = A( 7 , E) + 2^(7, E) + £ 2 A(- 7 , E), 
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/ U^r; E^ n )U^r; E')dr = 0, \E'\ > m, 
Jo 

POO 

/ U^r- E^ n )U 6 (r; E^)dr = 0, n ^ n' . 
Jo 

The normalization factor A n for the eigenf unctions of the discrete spectrum, 

POO 

A 2 n = / Ul(r;E (>n )dr, 
Jo 

can also be calculated (see, for example, [H]). It is interesting to note that we can explicitly 
verify the fulfillment of the relation A n = Q^\- For this purpose, we consider the integral 

/ C^(r; Et >n )V(r; E')dr = lim ^ ' g \ V = V{r; E), \E'\ < m, 

Jo E - 

where V(r; i?) is defined by (1701) . Using relations (ITT]) - (1T41) . we find 

f^U i {r-E^)V{r-E')dr = ^^ Wi . 
We now recall that V(r; ^, n ) — ^( r ; -^.n) and finally obtain that 

We also note that it follows from the inversion formulas that C € (E) = Q^ 2 {E). 
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